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ABSTRACT 
The  diabatic  heating  rate  at  lOOmb   is  determined  from  the  thermo- 
dynamic equation,  utilizing  previous  calculations  of  the  potential 
function  x  an<*  tne  "vertical  velocity"  oo,  for  the  period  1-5  April  1963, 
The  four  basic  terms  of  the  expanded  thermodynamic  equation;  storage, 
stream-velocity,  x~Poten^a^-   advection,  and  vertical  heat  transport  are 
examined  individually,  and  added  to  obtain  Q,  the  diabatic  heating  rate 
which  is  compared  to  the  results  of  previous  investigations.   The 
zonally-averaged  heating  rates  are  discussed  in  terms  of  stratosphere 
circulations  and  their  associated  physical  mechanisms. 
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1.    Introduction 

There  has  been  recent  interest  in  the  diabatic  heating  rates  and 
related  energy  changes  in  the  upper  troposphere  and  lower  stratosphere, 
notably  by  Davis  [1],  Julian  [2],  Martin  [4],  K.  Myakoda  [3],  and  others. 
Diabatic  heating  is,  of  course,  also  at  the  heart  of  the  problem  of 
large-scale  energetics.   However,  the  usual  approach  through  energetics 
yields  average  values  for  a  hemisphere  for  the  period  of  study,  whereas 
the  approach  of  this  study  takes  the  more  direct  route  of  obtaining 
12-hour  average  heating  rates  at  each  point  of  the  usual  octagonal 
(51  by  47)  analysis-prediction  grid  system  employed  by  the  Joint  Numer- 
ical Weather  Prediction  Unit,  and  other  groups  forecasting  large-scale 
weather  changes. 

Martin  [4]  demonstrated  the  feasibility  of  the  grid-point  approach 
by  obtaining  12-hour  mean  values  of  both  the  potential  function^  of 
the  wind  and  the  "vertical  velocity"  (*)=  -^.      Once  these  fields  are 
known  it  is  possible,  utilizing  the  thermodynamic  equation,  to  solve 
for  Q,  the  diabatic  heating  rate.  Martin's  model  is  described  in  his 
paper  [4]  and  only  those  aspects  required  for  this  paper  will  be  re- 
viewed here.   For  more  details  regarding  Martin's  model,  the  reader  is 
referred  to  his  paper. 

The  data  period  selected  for  the  calculations  was  1-5  April  1963. 
The  data  was  provided  by  the  Joint  Numerical  Weather  Prediction  Unit, 
Washington,  D.C.   This  particular  data  period  was  selected  in  an  effort 
to  reflect  the  relatively  large  heating  rates  and  the  associated  large 
meridional  flows  that  normally  occur  during  April,  particularly  in  the 
polar  stratosphere. 


2.    The  diabatic  equation 

If  one  expands  the  equation  of  thermodynamics,  the  result 


eat    jt 


p 


eap"       cP  (1) 

follows.   Here  Q  is  the  heating  rate  per  gram,  while  4>     and  X  give 
the  nondivergent  and  irrotational  parts  of  the  wind  in  accordance  with: 


y  =  kxv^+vX 


Here  the  fields  of  $>     and  X  were  obtained  preliminarily  by  Martin  [4] 
by  solution  of  the  Balance  Equation  and  the  X-2  equation  respectively 
(see  section  3).   When  the  equation  of  state  and  hydrostatic  equation 
are  employed  in  (1),  we  obtain  our  basic  working  equation: 


m^AtiY^Hi^i^-  * 


(2) 


The  fourth  term  of  (2)  may  be  transformed  into  the  more  convenient  form: 


*? 


n\l% 


"1   (3) 


where  (T*  =  -  -L  2&  but  is  also  the  negative  of  the  coefficient  of  U) 

P  \       /   \ 

in  the  final  form  of  (3).   Here  ( -rrr  I   and  f*r*l   are  centered  at  75mb 

V&PA*.    VApJi 

and  125mb  respectively,  that  is,  at  the  upper-half  and  lower-half  layer 
centers,  respectively.   The  term  -  (J"  of  (3)  has  a  positive  sign  since 
both  O-i  /tk)  and  the  contents  of  the  last  bracket  are  negative.   This 
agrees  with  the  intuitive  expectation  that  an  updraft  (w<o)  is  associated 
with  vertical-motion  cooling. 

The  factor  ^^/<^^/f\p  is  present  in  all  terms  on  the  left  side  of  (2)  , 


and  in  order  to  determine  its  value  at  lOOmb,  the  50mb  and  150mb  data 
levels  must  be  utilized.   Unfortunately  there  was  no  vertically  con- 
sistent 50mb  data  in  punched-card  form  immediately  available,  so  it  was 
necessary  to  derive  the  50mb  contour  data  completely  from  regression 
equations.   In  order  to  do  this  the  United  States  Navy  Weather  Research 
Facility  Regression  equations  of  D.  Lea  [5]  were  used.   These  equations 
have  the  following  form: 


z, 


=  K.+K^+l^ 


where  KQ,  K-^ ,  K2  represent  the  regression  coefficients.   It  should  be 
noted  that  the  regression  coefficients  KQ,  K-^,  K2  are  listed  by  ten 
degree  latitude  bands  from  ten  degrees  north  to  the  pole.   It  was 
necessary,  therefore,  to  devise  a  computer  program  converting  each 
Kn  ('(;)  to  a  value  of  K^  =  Kn  (i,j)  defined  at  each  grid  point  within 
the  51  by  47  octagon.   This  conversion  was  based  upon  a  subprogram 
giving  sin  (J)  in  terms  of  the  (i,j)  grid.   In  addition  Kn  (J))  was  made 
continuous  with  respect  to  latitude  by  averaging  the  values  across  the 
defining  latitude  boundaries  for  each  ten  degree  band,  while  retaining 
the  latitude-band  mean  value. 

Since  the  lOOmb  contour  analyses  during  the  period  of  study  were 
considered  to  be  geographically  consistent,  whereas,  those  of  tempera- 
ture were  not,  it  was  decided  again  to  obtain  T^qq  using  T-iqq  =  -JL.  ■  — 
in  (4).   With  this  modification,  the  regression  equation  is  expressible 
in  the  following  form: 

l-K,' 


where 


Kx  =   (K,  -  1) 


(&5Q  -  Z-^qq)   =   estimated  thickness  of  the  layer  100  to  50mb. 
and  the  term  ai*/^&r\p  appearing  in  (2)  has  been  expressed  in  the  form: 

£  fa  p  JU3  (6) 

In  (5)  and  (6) ,  (*cq  -  *ino^   ^s  estimated  from  the  modified  regression 
equation  (5)  ,  with  input  data  *ioO  anc*  (*100  ~  *150^  obtained  from 
copies  of  data  tapes  provided  by  the  Joint  Numerical  Weather  Prediction 
Unit.   The  parameters  £;loO  anc*  (*100  "  *150^  °^  (5)  an<^  (6)  were  used 
wherever  possible  in  order  to  utilize  as  much  vertically  consistent 
data  as  possible,  thereby  tending  to  minimize  errors  due  to  the  use  of 
the  regression  equations  (4)  or  (5).   All  computations  were  performed 
using  the  CDC  1604  computer. 

Equation  (3)  can  now  be  written  in  its  final  form: 


c, 


g^rvjQteU&tffrj  -j^'ifo^+kj  ^  x 


(7) 


15 


=  Q 


It  should  be  noted  that  2  =  D  +  £_,  where  *  is  the  standard  atmosphere 
height  and  D  =  £   -  *  is  the  height  anomaly  relative  to  that  of  the 
standard  atmosphere.   In  the  first  three  terms  of  (7)  %     is  independent 


of  the  variables  x,y,t  so  that  the  time  and  space  derivatives  become 


simply: 


1 


1 


V[(D^D^+(D^"D'^ 


etc.,  and  hence  £y,  can  be  omitted.   However  this  cannot  be  done  in  the 
vertical  motion  term  where  *  affects  the  vertical  stability. 
3.    The  irrotational  wind  and  vertical  motion  structure. 

As  pointed  out  in  section  1,  the  development  in  this  section  follows 
Martin  [4] ,  who  solves  diagnostically  for  x   for  his  x~2  equation  which 
is  written: 


(8) 


This  equation  is  solved  by  relaxation  using  an  iterative,  scan-procedure, 


subject  to  the  convergence  requirement: 

Then  with  x  assumed  to  be  known,  one  may  compute 


<  4.0 


cm 


v-v  o,  |v"X  =  -  (%0 


(9) 


Using  w  =  u)  (p)  ,  as  a  cubic  polynomial  in  p,  with  an  upper  boundary 
condition  to  =  0  at  p  =  0,  Martin  also  derives  an  integrated  form  of  (9) 
giving  w  =  or    '  at  each  grid  point. 

In  terms  of  finite  difference  notation  uA-1-00'  has  been  shown  [4] 

to  be:  (itfo)         , Z/W 

"J  .f   V  '        J  (10) 

Here   x  is  scaled  in  terms  of  centimeters  as  this  scaling  was  defined 


is  (8) 


From  (10),  (i),.,    '  takes  on  the  finite  difference  form: 


l> 


(11) 


with  x-j-j  assumed  zero  on  the  two  outermost  rows  of  the  octagon.   This 

was  also  the  boundary  condition  used  by  Martin  in  solving  (8) .   The 

47  by  51  octagonal  grid  mapping  program  used  in  the  computations  was  furnished 

on   tapes  by  the  National  Center  for  Atmospheric  Research.   Since  it  is 

not  possible  to  compute  >y  /\  on  the  "outer  boundary"  (Fig.  1),  it  was 

necessary  to  set  to  =  o  on  the  outer  boundary  throughout  the  computations. 
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Fig.  1 
The  boundary  conditions 
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Finite-differencing  grid 

The  normal  five-point  map  grid  1-2-3-4-0  was  used  one  row  within  the 
computational  boundary,  while  a  nine-point  grid  was  used  at  all  points 
further  inside.   The  nine-point  grid  simply  inserts  a  three-point  mean 
for  a  normal  difference;  for  example,  the  X3  "  Xi  difference  of  a  five- 
point  grid  is  replaced  with  the  weighted  difference: 

in  the  nine-point  grid.   The  actual  five-point  f inite-dif f erencing 


schemes  used  in  the  computations  are  as  follows  r?l9tfV€  to   tdfm> 


*V* 


l^4,  &iL 


n  v>^p]    2At 


vti 


VPX'V(sipK-^ 


Furthermore  the  resultant  w-jj  from  (11)  have  been  smoothed  twice  every- 
where interior  to  the  two  boundary  rows  just  mentioned,  giving  w-h  ,  the 
values  used  in  this  study. 
4.    Results 

With  the  fields  of  Xji  ar*d  u).  •  known,  as  well  as  \\jj-    and  observed 
values  of 


J. 


from  consecutive  12-hour  map-fields,  Q  has  been  determined  from  (7) 
The  first  term  of  (7) ,  called  the  "storage  term"  is 

I 


cjLfA- 


(12) 


The  second  and  third  terms  represents  the  12-hour  mean  horizontal  heat- 


transport  at  the  gridpoint  (i,j)  and  have  the  combined  form: 


h^K^s])+yM^H^ 


4 a3) 


'>  -  transport" 


X  -  transport 


The  fourth  term  in  (7)  is  the  vertical  heat  transport  ((^-transport) 


\F=>0    HOQ/  Vn6Q  "^sJ 


(14) 


The  four  terms  on  the  left  side  of  (7),  which  have  been  listed  in 
the  preceding  paragraph,  have  also  been  computed  individually  as 
functions  of  (i,j)  everywhere  within  the  computational  boundary.   Com- 
bined values  have  been  formed  to  give  resultant  octagonal  fields  of 
Q(i,j).   For  a  typical  field  of  Q  see  Fig.  3.   Zonal-averages  of  each 
of  the  five  fields  have  been  computed  by  programming  an  averaging 
routine  for  each  five  degree  latitude  band  (over  an  area  of  20-90 

degrees  latitude) . 

« 

Q,  and  the  preceding  four  terms  as  well,  have  units  of  gpm  which, 

for  later  comparison  purposes  must  be  changed  to  langleys.   Q  may  be 

determined  by  a  sum  of  terms  of  the  form: 


q  =  M 


<L-(3**) 


>0Q  Rd(U3^t 


(15) 


*  — 1        —  1—9 

Here  Q  is  in  ergs(gm)    (sec)  x  cm  ^  provided  all  terms  on  the  right 
side  of  (15)  are  in  c.g.s.  units,  AZ  =  %"^,  in  gp-cm  and  M,  nn  =  mass  of 


SOW 


Fig.  3.   Example  of  printout  of  Q  values  for  (06  GMT,  1  April  1963) 


the  column  150mb  to  50mb. 


^4 


(16) 


The  expression  within  the  braces  of  (16)  is  now  identical  to  (^TA^TL) 
in  units  of  (°K)(sec)~  .   To  get  a  corresponding  value  of  Q  per  gram 
per  day,  both  sides  of  (16)  are  multiplied  by  86,400  sec  (dayj,   and 
divided  by  102- the  latter  being  the  number  of  gms  in  lOOmb: 


'I 


QrrAQGtxAJ 


S$-\ZAl<i*l0-*  84.4QO 


9%c 


3py« 


sec 


Q  ■  6-298  x  [4.  <Azgpm]    •  •  •  • 


ly  (gm)  X  (day) 


The  units  of  ^  /<^t  [A&]  are  emphasized  as  gpm  (sec)-  because  this  is 
the  form  of  the  computer  output  units.   Thus  it  follows  f rom  ]<=  (AaSJ/^M 
that  a  geopotential  change  of  1  gpm  (gm)~  sec-1  is  equivalent  to  a 
heating  rate  of: 

6.298  ly  (gm)-1  (day)-1 

The  last  statement  is  typical  of  each  of  the  terms  in  (12),  (13),  (14). 
Zonal  averages  for  each  term  and  their  zonal  mean  sum  [Q]  of  (7)  have 
been  determined  and  are  listed  by  latitude  belts  in  Tables  1  through  5. 

# 

This  is  done  also,  for  purposes  of  comparison,  with  computations  of  [Q] 
made  by  Davis  [1],  who  used  a  climatological  or  standardized  atmosphere 
(Table  6).   Note  that  the  notation  [  ]  indicates  a  zonal  average  at  a 


10 


specified  latitude. 

In  Table  6,  the  nine  sets  of  results  obtained  here  are  averaged 
to  give  a  mean  daily  average  of  the  storage,  the  horizontal  transport 

terms  of  both  the  stream  and  potential  types,  the  vertical-heat 

•  _  I     _  -i 

transport  and  the  resultant  Q,  each  in  ly(gm)    day.    Davis1  results 

have  been  centered  at  117. 5mb,  but  presumably  this  slight  shift  of 

levels  is  insignificant.   In  regard  to  the  units  of  ly(mb)--1-  day 

attributed  to  Davis,  it  should  be  recalled  that 


l^jL  =    Lsgj 


whereas  out  unit  of  heat  is: 


liiL.  =   I  cat 

9™         %0  iy-viiS 


Thus  the  two  heating  rates,  of  this  paper  and  of  Davis,  are  very  nearly 
compatible.   The  essential  results  are  presented  here  as  meridional 

cross-sections  contained  in  Tables  1  through  6. 

4 
In  Table  6b,  Davis'  computations  for  mean  daily  Q  are  based  upon 

a  climatological  atmosphere  at  rest.   It  will  be  seen  that  Q  should 

depend  upon  the  nature  of  the  meridional  motion  pattern,  especially 

the  field  of  vertical  motion,  in  addition  to  any  static  distribution 

of  radiating  material.   The  term  "Balance"  is  taken  from  Davis'  paper 

and  refers  to  balance  requirements,  on  the  basis  that 

-<£.r)=  q  =  tpjVf  ^T+v^T-cruo] 

with  the  "storage"  very  small  (by  several  orders  of  magnitude)  both  in 
the  present  study  and  in  Davis'. 
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5.    Conclusions:   stratospheric  circulations  and  associated  physical 
mechanisms 

As  already  noted,  the  zonally-averaged  storage  term  is  very  small 
In  addition  the  stream-  and  potential-advections ,  averaged  over  five- 
degree  latitude  bands  are  smaller  by  approximately  two  orders  of  magni- 
tude than  [Q] .   Thus  to  a  good  approximation,  valid  everywhere  at 
least  as  to  sign,        ("  #  ")        C  "1 

and  is  of  the  proper  order  of  magnitude  as  well. 
By  the  equation: 


(17) 


it  follows  that  [w]oq]  >0  occurs  when  [^V-^]  rePresents  convergence. 
Now  convergence  occurs  near  the  pole,  since  the  expression  for  the 
zonally-averaged  divergence 


Lv'%d    U  a)  UwJ^T-*"  (18) 


—■■in  * 


with  the  second  term  (which  represents  the  convergence  of  meridians) 
becoming  dominant  [Martin,  4]  at  latitudes  lD  >  75N.   At  lower  latitudes 

LV'Voo>  [^ao>  (^/^] 

is  a  valid  approximation;  that  is,  the  spherical  divergence  is  only 
slightly  affected  by  convergence  of  meridians.   Hence,  where  Ico-,  ^Jj  is 
negative,  there  is  divergence,;   ifh,QQlis  positive,  there  is  convergence, 
Consequently  at  lOOmb,  we  have  centers  of  convergence  at  the  pole  and 
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latitude  35,  and  of  divergence  at  latitude  65. 

In  inferring  a  stratospheric  circulation,  one  may  start  with  an 
interhemispheric  meridional  circulation  proposed  by  Kellogg  [6], 
especially  applicable,  in  some  of  its  aspects,  to  the  winter-summer 
reversal  in  the  stratosphere.   The  octagonal  grid,  together  with  the 
associated  boundary  conditions,  is  not  able  to  resolve  the  effect  of 
the  cross-equatorial  flow  to  at  least  latitude  20N  latitude.   However, 
because  of  the  decrease  of  mass  in  the  winter  stratosphere,  cross- 
equatorial  flow  must  occur  near  the  time  of  the  vernal  equinox.   This 
evidently  is  responsible  for  the  convergence  center  near  latitude  35N. 

Likewise  with  the  return  of  the  sun  to  the  Northern  Hemisphere, 
ozone  regeneration  goes  on  rapidly,  followed  by  ozone-heating, 
especially  at  45  km.   According  to  Pressman  [7],  an  ozone  radiative 
heating  maximum  proceeds  progressively  from  30N  latitude  to  the  pole 
between  the  dates  15  March  and  15  May.   For  time  of  this  study,  1-5 
April,  1963,  one  might  expect  maximum  ozone  heating  at  45  km  at  or  near 
60N  latitude.   This  effect  is  so  marked  as  to  cause  a  strong  forced 
outdraft  (V'^V^O)  at  60N  latitude,  and  at  an  elevation  of  45  km.   In 
order  that  three-dimensional  mass  divergence  be  within  reasonable 
bounds,  there  must  be  a  sympathetic  updraft  from  below.   This  is  the  [to] 
measured  at  the  lOOmb  level  and  presumably  exists  at  levels  between 
lOOmb  and  45km.   Also  according  to  Martin's  model,  sizeable  co-values  of 
the  same  sign  still  persist  at  the  200mb  level,  but  are  smaller  than 
toiOO  in  units  of  cm  sec--'-.   The  level  of  "return  flow"  in  these  cells 
cannot  be  designated  from  presently  known  studies. 
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Fig.  4.   Proposed  meridional  flow  patterns 

early  in  the  spring-reversal  period 


The  "maximum"  convergence  in  the  vertical  and  associated  downdrafts  at 
the  pole  at  this  particular  time  of  ozone  activity  should  then  be 
viewed  as  a  "spill-over"  effect  from  65N  latitude.   The  maximum  con- 
vergence at  or  near  latitude  35  North,  which  mirrors  a  zonally-averaged 
100-mb  maximum  values  of (a)  at  this  same  latitude  is  postulated  as  the 
primary  result  of  cross-equatorial  flow.   Considering  the  level  at  which 
the  greatest  polar-vortex  intensification  in  the  Southern  Hemisphere 
occurs,  the  level  of  the  35N  latitude  maximum  of  convergence  has  also 
been  placed  at  45km. 

Fig.  4  together  with  the  above  discussion  "explains"  the  mean 
meridional  flow  pattern  in  the  stratosphere  during  the  period  of  this 
study.   One  may  offer  as  conjecture  that  as  the  stratospheric  summer 
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advances  (say  by  June  15) ,  the  maximum  convergence  will  have  advanced 
to  45N  latitude  while  the  zone  of  maximum  divergence  associated  with 
most  effective  ozone-heating  near  45km  advances  to  the  pole,  there  to 
become  associated  with  divergence  at  higher  elevations. 

As  contrasted  with  the  Q  of  this  study,  Davis'  results  show  only 
moderate  cooling  at  every  latitude  due  to  predominance  of  long-wave 
radiation  over  heat  sources.   Two  features  by  which  April,  1963 
appears  to  differ  from  climatology  are  shown  in  Figs.  5  and  6,  both 
taken  from  analyses  of  the  Free  University  of  Berlin  [8].   The  latter 
figure  shows  a  pre-existing  zone  of  vortices  at  latitude  60  North, 
as  of  1  April,  1963,  1200  GMT.   The  other  figure  mentioned  above  shows 
that,  as  compared  with  a  completely  different  April  (1962),  pressure 
was  higher  over  the  pole  but  lower  than  in  1962  in  a  vortex  ring  at 
60N  latitude.   Fig.  6  shows  a  possible  reason  for  this:   strong  baro- 
clinicity  on  the  equatorward  side  of  the  vortices,  with  'zonal 
generation' 


Q  T 


/ 


computed  over  45-70  degrees  North  latitude  favoring  the  continuation  of 
such  a  baroclinic  zone  and  westerly  winds.   In  fact,  according  to  the 
analyses  of  the  Free  University  of  Berlin,  an  easterly  regime  with  an 
associated  polar  anticyclone  did  not  occur  until  May  15,  1963. 

Finally,  how  may  one  reconcile  the  zonally-averaged  equation: 


Q-(-c,xrui)^o 


(19) 


so  that  both  horizontal  flows  and  storage  are  relatively  unimportant 
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Fig.  5.   lOOmb  height  changes,  April  1963 
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as  determined  here.   Equation  (19)  implies  that  with  —  C*pU-)J>0   (d 
draft  motions) , IqJ<0  is  the  result.   Such  a  result  is  reasonable  on 
physical  grounds  by  noting  that  descending  air  is  warmed  dry  adiabati- 
cally,  and  in  addition  one  may  expect  some  increased  diabatic  solar 
heating  acting  upon  an  ozone  increase  in  lower  stratospheric  elevations 
However,  one  immediate  opposing  result  is  that  the  CO2  in  the  warmed 
air  will  emit  at  a  greater  rate,  tending  to  re-establish  the  heat 
balance.   The  same  argument  applies  in  reverse  for  updraft  motions. 
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